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ABSTRACT 


A  combination  of  Interactive  boundary  layer  and  stability  theories  has  been 
used  to  Investigate  the  reasons  for  the  Instability  of  laminar  separation 
bubbles  on  the  leading  edge  of  thin  airfoils.  It  Is  shown  that  transition 
plays  an  Important  role  and  is  likely  to  preclude  the  existence  of  long 
separation  bubbles  and  their  supposed  Instability. 
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1.0  INTRODUCTION 


The  laminar  boundary  layer  on  an  airfoil  grows  from  the  stagnation  point  with 
a  favorable  pressure  gradient  which  causes  the  flow  to  accelerate  and  Is  then 
subjected  to  an  adverse  pressure  gradient  which  can  cause  separation  with  sub¬ 
sequent  reattachment.  The  resulting  bubbles  are  common  on  thin  airfoils  where 
the  adverse  pressure  gradient  can  be  sufficiently  strong  to  cause  the  flow  to 
separate  even  at  small  angles  of  attack  and  are  Important  because  of  their 
association  with  the  phenomenon  of  stall.  The  nature  of  the  phenomenon  Is 
known  to  depend  on  the  Reynolds  number  based  on  the  radius  of  the  leading  edge, 
as  will  be  shown,  the  length  of  the  separation  bubble  can  grow  to  Influence 
the  location  of  transition  and,  on  occasions,  transition  can  occur  within  the 
bubble.  The  angle  of  attack  Is  also  known  to  be  Important  and  can  cause  the 
separated  region  to  grow  until  at  some  angle  of  attack  the  bubble  bursts  and 
stall  sets  In.  This  sequence  of  events  from  the  first  appearance  of  separa¬ 
tion  on  the  upper  surface  to  stall  Is  complex  and  requires  clearer  understand¬ 
ing  than  Is  presently  available. 

The  prediction  of  separation  bubbles  on  airfoils  has  been  studied  by  a  number 
of  Investigators.  An  Important  contribution  was  made  by  Briley  and  McDonald 
[1]  who  used  an  Interactive  boundary-layer  approach  and  solutions  of  the 
Navler-Stokes  equations  to  examine  the  case  of  a  comparatively  thick  airfoil 
where  the  separation  region  occurred  around  midchord  and  was  approximately  10% 
chord  In  extent.  Crlml  and  Reeves  [2]),  Kwon  and  Pletcher  [3],  Cebecl  and 
Schlmke  [4]  and  Carter  and  Vatsa  [5]  have  tackled  essentially  the  same  problem 
with  Interactive  boundary-layer  theory.  In  all  cases  the  location  of  transl- 
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tlon  was  either  assumed  to  correspond  to  laminar  separation  or  was  computed  by 
an  empirical  formula.  The  work  of  Cebecl  and  Schlmke  also  examined  the  Influ¬ 
ence  of  the  location  of  transition  and  showed  that  reattachment  and  transition 
were  related.  Attempts  to  perform  calculations  with  the  experimentally 
reported  transition  location,  which  occurred  further  downstream  than  those 
considered  above,  revealed  a  tendency  for  the  reattachment  location  to  move 
rapidly  downstream  with  the  number  of  sweeps  used  In  the  Interactive  proced¬ 
ure.  This  apparent  Instability  of  the  separation  bubble  Is  examined  further 
In  this  report  together  with  Its  relationship  to  the  location  of  transition. 


The  present  approach  can  be  described  In  two  parts.  First  we  will  make  use  of 
linear  stability  theory  and  the  en-method  to  predict  transition  based  on 
calculated  velocity  p  oflles.  The  validity  of  this  procedure  has  been  demon¬ 
strated  by  Cebecl  and  Egan  [6]  and  calculations  are  presented  to  confirm  that 
It  Is  appropriate  for  the  particular  flows  under  Investigation  here.  Secondly, 
we  follow  the  approach  of  Cebecl  et  al.  [7]  and  examine  the  leading-edge  sepa¬ 
ration  bubble  on  a  thin  airfoil  as  a  function  of  angle  of  attack  and  for  a 

5 

Reynolds  number  of  10  .  This  systematic  study  has  been  arranged  to  allow  us 
to  examine  carefully  the  relationship  between  the  growth  of  the  separated 
region  and  transition. 


2.0  CALCULATION  OF  TRANSITION 


In  a  recent  study,  Cebecl  and  Egan  [6]  performed  calculations  of  steady  flows 
over  and  downstream  of  bumps  Identical  to  those  examined  experimentally  by 
Fage  [8].  The  shape  of  the  bump  was  represented  In  the  calculations  by 


1  -  12  (jp  -  16  (jp 
1  -  12  (|)2  16  (|)3 


<  x  <  0 


0  <  x  <  | 


where  h  and  B  denote  the  height  and  width  as  shown  In  Figure  1.  The  calcula¬ 
tion  method  was  based  on  that  of  Cebecl  et  al.  [7]  and  solved  the  boundary- 

layer  equations  In  an  Inverse  mode  with  successive  sweeps  over  the  body.  The 

o 

edge  boundary  condition  was  written  as  the  sum  of  the  Invlscld  velocity  ug(x) 

and  a  perturbation  velocity  6u  (x),  that  Is, 

e 


y  =  6, 


ue(x)  =  ue(x)  *  6ue(x) 


and  4ue(x)  was  obtained  from  the  Hilbert  Integral  given  by 

x 

/  \  1  ,  b  d  ,  do  ,  ^ 

u  (x)  =  -  J  ~r~  ( U  6*)  -  (3 

e  '  w  da  '  e  '  x  -  a 

a 

with  the  Interaction  region  confined  between  x  and  xu . 

a  b 

o 

fhe  above  Interactive  boundary-layer  procedure  with  up(x)  =  1  was  used  to 
compute  the  boundary-layer  characteristics  Including  the  velocity  profiles  and 

tl 

wall  shear  stress  parameter  f  defined  by 

w 


for  the  conditions  Investigated  by  Fage.  A  sample  of  the  results.  In  terms  of 

"  5 

f  ,  are  Included  In  Figure  1  for  a  Reynolds  number  of  4.375  x  10  per  foot 

and  for  three  bump  heights.  Here  the  Reynolds  number  per  foot  Is  defined  In 
terms  of  the  measured  freestream  velocity  u^  at  the  position  of  the  center¬ 
line  of  the  bump  but  for  the  undistorted  surface.  The  figure  shows  that  the 
wall  shear  parameter  decreases  Immediately  prior  to  the  bump,  rises  rapidly 
with  the  favorable  pressure  gradient  Imposed  by  the  upstream  surface  of  the 


(comp)  =  5.29 


=  0.0310" 


x.  (exp)  =4.58 
tr 


(comp) = 3.85  » 


=  0.0525" 


xtr  (exp)  =  2.92 


(comp)  =  2.95 


=  0.0620" 


ngure  i.  variation  or  wan  snear  parameter 
for  constant  Reynolds  number. 

bump,  reaches  a  maximum  and  decays  rapidly  t 
Ing  as  the  Influence  of  the  bump  diminishes. 
Is  to  Increase  the  magnitude  of  the  maxima  a 
tlon  with  corresponding  Increase  In  Its  grad 


Figure  1  and  Table  1  also  show  measured  and  calculated  locations  of  transition 
with  the  latter  obtained  from  the  e  -method  and  the  calculated  velocity  pro¬ 
files.  This  method  stems  from  the  work  of  Smith  and  Gamberonl  [9]  and  Van 
Ingen  [10]  and  Is  based  on  linear  stability  theory.  It  assumes  that  transi¬ 
tion  starts  when  a  small  disturbance  Is  Introduced  at  a  critical  Reynolds  num¬ 
ber  and  Is  amplified  by  a  factor  of  eP.  For  given  velocity  profiles,  the 
Orr-Sommerfeld  equation  is  solved  and  stability  properties  are  examined.  The 
amplification  rates  (-a^)  are  computed  as  a  function  of  x  for  a  range  of 

discrete  values  of  the  frequency  u>  and  transition  Is  assumed  to  occur  when 
-fa  dx  n 

e  1  reaches  a  value  equal  to  e  where  n  Is  around  9.  In  the  present  case 
the  profiles  were  available  from  the  Interactive  boundary-layer  calculations 
and  the  same  version  of  the  Box  scheme  was  used  to  solve  the  stability  equa¬ 
tion  with  a  continuation  method  to  obtain  the  eigenvalues  In  regions  of  rap- 

II 

Idly  changing  f  and  In  regions  of  separated  flow. 

w 

The  agreement  between  measured  and  calculated  transition  locations  Is  shown  on 
Figure  1  and  Table  1  to  be  within  experimental  uncertainty  and  similar  results 
were  reported  by  Cebecl  and  Egan  for  the  much  wider  range  of  configurations 
and  Reynolds  numbers  Investigated  by  Fage.  It  Is  clear  that  the  location  of 
transition  moves  upstream  with  Increasing  bump  height  and  that  the  length  of 
the  separated  region  Increases.  These  two  characteristics  are  also  to  be  found 

II 

In  the  f  distributions  associated  with  the  leading  edge  region  of  thin 

w 

airfoils  as  discussed  below. 

II 

Figure  2  shows  the  fy  distribution  and  the  corresponding  external  velocity 
distribution  for  the  leading  edge  of  a  thin  airfoil.  The  result  Is  similar  to 

Table  1.  Influence  of  Freestream  Velocity  u^c  on  Transition, 

h  =  0.0620" 

Transition  Length,  L  (ft) 
measured  *~calc 


u  (ft/sec) 
Ic _ 


3.0  SEPARATION  BUBBLES  ON  THIN  AIRFOILS 


following  the  approach  of  Cebecl  et  al.  [7],  we  consider  a  thin  ellipse  In 
nose-centered  coordinates 

(4) 

a  b 

at  an  angle  of  attack  of  a  and  In  a  uniform  stream  of  velocity  u  .  Attention 

00 

Is  directed  to  the  nose  region  where  the  ellipse  Is  approximated  locally  by  a 
"nose-fitting"  parabola,  and  an  expression  Is  derived  for  the  external  velocity 
distribution  ue. 


ue  =  u  (1  +  t) 


5  *  l. 


'1  ♦  V 


Here  t  denotes  the  thickness  ratio  b/a,  E  corresponds  to  a  reduced  angle  of 

o 

attack,  a/t,  and  the  parameter  E  Is  a  dimensionless  distance  from  the  nose 

2  2 

related  to  the  x-  and  y-coordlnates  of  the  ellipse  by  x  +  a  =  1/2  at  E  .  y  = 

2 

at  E-  The  parameter  E  Is  related  to  the  surface  distance  s  by 

s  =  at2  (1  ♦  E2)1/2dE  (6) 

0 

The  Investigation  of  Cebecl  et  al.  [7]  made  use  of  the  external  velocity  dis¬ 
tribution  given  by  Eq.  (5)  and  showed  that  the  laminar  boundary  layer  near  the 
leading  edge  was  well  behaved  and  unseparated  If  EQ  <  E$  =  1.16,  although 

there  was  significant  adverse  pressure  gradient.  At  higher  values  of  E  . 

o 

however,  separation  occurred  with  an  associated  singularity  and  required  the 
use  of  an  Interactive  theory  to  link  the  viscous  and  Invlscld  flows.  With 

this  theory,  solutions  were  obtained  for  separation  bubbles  at  R(E  2  u  a/v)  = 

6  00 
2  x  10  and  for  t  =  0.1  but  reattachment  occurred  In  a  very  limited  range  of 

the  reduced  angle  of  attack.  For  EQ  >  1.218,  calculations  broke  down 

shortly  after  the  flow  reversal  In  the  boundary  layer  and  the  subsequent 

studies  of  Stewartson,  Smith  and  Kaups  [11]  led  them  to  suggest  that  a  dramatic 

switch  to  another  separated  form  of  motion  can  occur. 


Similar  calculations  have  been  performed  for  a  Reynolds  number  of  10  ,  for 
which  we  can  be  sure  that  the  transition  occurs  In  the  region  downstream  of 


the  separated  flow,  and  are  presented  In  Figures  3  and  4  with  Increase  of  the 

reduced  angle  of  attack  Implying  an  Increase  In  the  strength  of  the 

adverse  pressure  gradient.  Consistent  with  the  observations  of  Cebecl  et  al. 

[7],  the  region  of  separated  flow  was  found  to  Increase  In  extent  with  l  and 

o 

at  l  =  1.296  revealed  a  tendency  to  expand  slowly  with  each  sweep  (Fig.  4) 
o 

with  the  tendency  becoming  bigger  at  l  =  1.298  (Fig.  5).  This  Instability  of 

o 

the  bubble  may  have  a  counterpart  In  the  bump  flows  of  Section  2  since,  with  a 
very  severe  adverse  pressure  gradient  corresponding  to  a  larger  bump  than  those 
considered  by  Fage,  It  may  be  expected  that  laminar  flow  will  undergo  transi¬ 
tion  before  reattachment.  This  possibility  can  be  tested  with  the  help  of  the 
stability  theory  described  In  the  previous  section. 


The  application  of  the  en-method  to  the  leading-edqe  flows  of  Figure  3  led 

tc  transition  locations  Identified  In  the  figure.  They  exhibit  the  same  trend 

noted  In  connection  with  Figure  1  In  that  transition  moves  forward  with 

Increasing  reduced  angle  |  which  Is  analogous  to  bump  height  h.  It  Is  of 

particular  note  that,  as  £  tends  to  the  value  of  1.296  for  which  Insta- 

o 

blllty  has  been  observed,  the  transition  location  moves  towards  and  Inside  the 
separation  bubble  at  j-  =  1.292.  With  further  Increase  In  the  reduced 
angle  of  attack,  the  transition  location  moves  further  Inside  the  separation 
bubble.  At  =  1.296  (see  Fig.  4),  the  transition  location  moves  upstream 
with  each  sweep.  These  results  Imply  that  the  real  flow  will  become  turbulent 
and  have  a  shorter  recirculation  region  which  Is  consistent  with  experiments. 
It  also  suggests  that  there  Is  little  merit  In  expending  effort  to  calculate 
the  large  laminar  separation  bubbles  which  would  be  obtained  with  larger 
reduced  angles.  This  observation  Is  likely  to  be  Independent  of  the  use  of 
Interactive  boundary  layer  or  Navler-Stokes  procedures. 
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